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Seiberg and Witten's proposed solution of N = 2 SQCD with N c = 2 and Np = 4 is 
known to conflict with instanton calculations in three distinct ways. Here we show how 
to resolve all three discrepancies, simply by reparametrizing the elliptic curve in terms of 
quantities and u rather than r and u = (Tr A 2 ). SL(2,Z) invariance of the curve is 
preserved. However, there is now an infinite ambiguity in the relation between and r 
and between u and u, corresponding to an infinite number of unknown coefficients in the 
instanton expansion. Thus the reinterpreted curve (unlike the cases Np < 4) no longer 
determines the quantum modulus u as a function of the classical VEV a. 
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1. Review of discrepancies between the Seiberg-Witten solutions and instanton 
calculations 



The solutions of N = 2 supersymmetric QCD with gauge group SU (2) proposed by 
Seiberg and Witten |JJ have recently been tested against first-principles instanton calcula- 
tions, up to the two-instanton level @-ff]. As reviewed below, these tests have resulted in 
perfect agreement when Np, the number of quark hypermultiplet flavors, is < 2. On the 
other hand, interesting discrepancies have emerged for the cases Np = 30 and Np = 4 
0. In this paper, we present our conjecture for how to resolve these discrepancies for 
Np = 4. Just like the conjectured resolution for Np = 3 0, we will not actually need to 
modify the elliptic curve that governs the solution (Eq. (16.35) of ) . Instead, we will 
change the interpretation of the parameters r and u that appear in this expression; they 
will have a different (and highly underdetermined) relation to physical observables of the 
microscopic theory than that proposed by Seiberg and Witten. In particular, in sharp 
contrast to the cases Np < 4, the reinterpreted curve for Np = 4 no longer determines the 
quantum modulus 

u = (TrA 2 ) , (1) 

as a function of the classical VEV a of the adjoint Higgs A = A a r a /2. 

The solutions given in make precise predictions for all multi-instanton contribu- 
tions to the low-energy physics. The functional form of these contributions is tightly 
constrained by holomorphy, the U{1)r anomaly, renormalization group (RG) invariance, 
and the Matone relation |6|-|8[] between the prepotential jF(a) and the modulus u(a). For 
Np > 0, additional constraints come from flavor symmetries, the decoupling limit for 
heavy flavors, and a discrete 2Z<i symmetry forbidding odd-instanton contributions when 
any quark mass vanishes. As we showed in Refs. [p|,p|,P,pX]| , all these constraints are built 
into the instanton calculus as well.0 Therefore, in testing the Seiberg-Witten solutions 
against explicit instanton calculations, one can "mod out" these constraints, and focus on 
a basic subset of nontrivial numerical predictions extracted from the elliptic curves at each 



1 The all-instanton-orders proof of the Matone relation given in fl~0[] was ostensibly limited to 
the case Np = 0. However, the reader can easily verify that the proof goes through for arbitrary 
numbers of (massless or massive) bypermultiplets. The key point here is that Eqs. (20)-(21) of 
[1C] can be directly extended to Nf > 0, see Eq. (7.20) of ||. The Matone relation can also be 
understood as a Ward identity for (broken) superconformal invariance [11]. 
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order in the instanton expansion. These numbers are then compared against the results of 
a finite-dimensional integration over the (multi-)instanton super-moduli. 



N F 



parameter 
fixing 



first nontrivial prediction first nontrivial prediction 
of curve (> 1 massless quark) of curve (all m q > ) 



1-loop 

1 1-loop 

2 1-loop 



1-inst 



N/A 



2-inst 



2-inst 



2-inst 



2-inst 



3 l-loop+2-inst 

[Eq. (4)] 

4 1-loop 

+ all-even-insts 



none if all m q = 0; 
otherwise 4-inst 



4-inst 



3-inst 



3-inst 



[Eqs. (6), (42)] 



Table 1. For each Np < 4, the second column gives the order in the semiclassical 
expansion at which the input parameters and u (or and u for Np = 
4) in the Seiberg-Witten curves are fixed in terms of physical quantities in the 
microscopic theory. The third and fourth columns (respectively, at least one 
massless quark, and all quarks massive) then give the order in the semiclassical 
expansion at which a nontrivial numerical prediction is first made by the curve, 
which can be compared against an instanton calculation ("nontrivial" is defined 
in the text) . To date, no numerical tests beyond the 2-instanton level have been 
performed, so that the models with N F = 0, 1, 2 are on a firmer footing than the 
ones with Np = 3, 4. 

In formulating these tests, and evaluating their outcome, it is crucial to distinguish 
between what is input (i.e., the fixing of conventions) versus what is output (i.e., a definite 
numerical prediction). Since this distinction will be key to resolving the discrepancies for 
Np = 4, we should first review the models with Np < 3 ; see Table 1 for a summary. The 
input parameters for the elliptic curves with < Np < 3 are the dynamically generated 
scale A^, and the modulus u, as well as the Np quark masses. (Notation: henceforth 
we will denote by u the parameter in the Seiberg-Witten curves, and reserve the symbol 
u always to mean (Tr A 2 ).) The curve defines an implicit choice of regular ization scheme 
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for Ajy^. For any physical quantity the n-instanton contributions are proportional to 
(A^) n ( 4-ArF ) ; in order to test the specific numerical predictions of the curves against 
instantons, one needs to know how to relate to (say) the analogous quantity A^ in 
the Pauli-Villars (PV) scheme which is most natural for instanton calculations [12|. As 



explained by Finnell and Pouliot []T3|] and reviewed in Sec. 2 below, this "dictionary" is 
fixed at the 1-loop level, and reads: 

(AZ) 4 ~ Nf = WJ 4 -^, JV F = 0,1,2,3. (2) 

Similarly, for < Np < 2, a discrete symmetry on the u plane combined with a semiclas- 
sical analysis of the singularities gives simply ]]J 

u = u , N F = 0,1,2 . (3) 

Equations (H)-(H) are the input; the predictions for J-(a) and u(a) first lie at the 1-instanton 
level for Np = 0, and at the 2-instanton level for Np = 1,2 due to the %2 symmetry 



mentioned above." The tests for Np = have been performed at the 1-instanton [13 
and at the 2-instanton || level, and the tests for Np = 1,2 have been carried out at the 
2-instanton level @-|4|], all with exact agreement. 

Next we consider the case Np = 3. Here Seiberg and Witten (in Sec. 14.1 of 0) 
continue to assume the equality (||), even though it is no longer required by a symmetry 
argument. Indeed, for Np = 3, the Z^-Np symmetry that acts on the complex u plane 
is trivial 0; consequently Eq. (||) can be generalized to u = u — uq ■ (A| w ) 2 where uq is a 
numerical constant. (This possibility is noted by Seiberg and Witten, but not exploited.) 
A nonzero value for uq has been extracted from a 2-instanton calculation by Aoyama et 
al; they find § 

u = u-uo-iAD 2 , w o = -^, N F = 3 (4) 

contradicting the naive assumption (|3|). It is natural to hypothesize that the solution 
to the Np = 3 model is still given by the Seiberg- Witten curve (Eq. (14.5) of 0), with 
the substitution (0J) rather than (|3|). In sum, for Np = 3, both the 1-loop level and the 



2 When quark masses are nonzero the models with Nf =1,2 have 1-instanton contributions 
too; however, due to built-in heavy-quark decoupling properties, these contributions are pegged 
to the 1-instanton term for Nf = and so are not independent tests 0,0]. 
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2-instanton level should be considered as input due to the extra degree of freedom m ; the 
first testable predictions extracted from the curve then lie at the 3-instanton level (or, if 
one of the quarks is massless, at the 4-instanton level). 

In this paper we will focus on the particularly interesting case Np = 4. Here, unlike 
the previous cases, the proposed solution and the instanton calculation disagree at the level 
of the low-energy effective U(l) Lagrangian £ e ff @- In fact, this disagreement can already 
be seen in the massless model. In this case both the /3-function and the U(1)r anomaly 
vanish identically, so that the microscopic coupling g± and ^-parameter 64 (which can no 
longer be rotated away) assemble into a single scale-independent holomorphic coupling 
t = + ^ • Furthermore is simply the Lagrangian of a classical N = 2 free field 
theory; its overall normalization, t^, enters the BPS formulae for the dyon masses 
(The superscript (0) will denote the massless case.) In Section 10 of Seiberg and Witten 
make the strong additional assumption that, thanks to the absence of a running coupling 
constant, the effective U(l) coupling equals the classical 577(2) coupling, 

t<% = r, (5) 

which implies no quantum (perturbative or non-perturbative) corrections to r. Instead, as 
we demonstrated in Ref . || , a first- principles instanton calculation gives! 

r W = i^c°)"( a ) = T + 1 c n q n , q = exp(zTrr) . (6) 

71 n=0,2,4... 

In particular a nonzero 2-instanton contribution C2 = — 2^35 was calculated in Ref. ||. 
The 1-loop perturbative constant Co, while not considered in Ref. ||, will turn out to be 
crucial to our proposed resolution. We calculate below that Co = 4 log 2 in the PV scheme.^ 
A second disagreement in the Np = 4 model arises when at least one of the four 
hypermultiplets has a nonzero mass, say 7^ 0. In the double scaling limit defined by 
—* 00 together with — ► in a specific way reviewed below, the heavy flavor decouples, 
and the model is supposed to flow to the iVi? = 3 theory. This requires the identification 



3 Note that we use q for the 1-instanton factor rather than for the 2-instanton factor as in Q. 

4 If the classical exactness assumption (|5j) were correct, there would be no need to introduce 
a regularization scheme, since scheme dependence is a one-loop effect. As discussed in Sec. 3, 
although the Nf = 4 model is a finite theory, the microscopic SU(2) coupling 54 must still be 
defined in a particular scheme; that is because the finiteness is due to cancellations between 
individually UV divergent graphs. 
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of the dynamical scale A™ for Np = 3, with the parameters of the Np = 4 theory. As 
explained below, working in the PV scheme and using the dictionary (^), one obtains: 

A|w = 4A pv = 4TO4 exp(-87r 2 /(^ v ) 2 ) , (7) 

where g\ w is the microscopic PV coupling in the 4-flavor model. In contrast, the relation 
given in [1]], 

Ai w = 64m 4 exp(-87r 2 /(d w ) 2 ) , (8) 

involves a proportionality constant of 64 rather than 4. Again, if the quantum corrections 
to Eq. (|D were absent, there could be no distinction between the PV coupling constant 
g\ Y and the "classical" g™ of Seiberg and Witten, hence no accounting for the factor of 
16 mismatch between Eqs. (Q) and (^). 

Thirdly, even with the identification (|]), the specific Np = 4 solution proposed in 
flows to the uncorrected version of the Np = 3 model which fails to incorporate the shifted 
definition of u given in Eq. (^). 

Our principal aim is to explain how to resolve all three of these Np = 4 discrepancies 
in a simple way, through a reinterpretation of the quantities r and u that enter into the 
massive Seiberg- Witten curve. In particular, rather than being modular forms of the 
microscopic £77(2) parameter r, the coefficients of the massive curve will be functions 
of the effective massless U(l) coupling defined by the all-even-instanton series (|6|). 
(Obviously the factor of 16 between Eqs. (0) and (H) will be automatically accounted for 
by the exponentiation of Co = 4 log 2.) This redefinition preserves the important SL(2, ZZ) 
invariance of the elliptic curve, as well as the stringent residue condition described in Sec. 17 
of . However, it also introduces an infinite ambiguity into the solution, parametrized by 
the infinite number of as- yet- undetermined numerical coefficients C4, cq, . . . , in Eq. ([|) (a 
similar series, Eq. (42), relates u and u). These issues are discussed in Sec. 3 below, which 
is devoted to the case Np = 4, and also contains our conclusions. But first, in Sec. 2, we 
lay some necessary groundwork in the cases Np < 4. In particular we review Weinberg's 
matching prescription between high- and low-energy gauge theories |14]| , as this formalism 



lies at the heart of the physics. In the process we will also specify the perturbative, one- 
loop contributions to the prepotential for Np < 4, correcting some incomplete expressions 
in the literature. 

The main result of this paper is the reinterpreted Np = 4 massive curve which agrees 
with all available perturbative and multi-instanton calculations. Although the input pa- 
rameters of this curve receive contributions from all even numbers of instantons, the curve 
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does contain definite numerical predictions at the 3-instanton level (assuming nonvanish- 
ing quark masses) that can, in principle, be tested against a semiclassical calculation. We 
anticipate that similar reinterpretations need to be made in the general class of N = 2 
models with gauge group SU(N C ) and Np = 2N C for which the /3-function vanishes. (In- 
deed, discrepancies at the 1-instanton level have been claimed for N c = 3, Np = 4, 6 by 
Ito and Sasakura fl5fl .) 



2. Weinberg's matching prescription and the Seiberg-Witten regularization 
scheme for Np < 4 

The physics of N = 2 SQCD utilizes — in two distinct but equally important ways — a 
matching prescription between a "high-energy" and a "low-energy" gauge theory. On the 
one hand, the SU (2) gauge group spontaneously breaks down to U{1) as the adjoint Higgs 
A acquires a complex VEV ( A) = ar 3 . For energy scales E <C Mw = v8|a|, the dynamics 
is governed by a nonrenormalizable Wilsonian effective action with U(l) gauge invariance, 
formally obtained by integrating out the heavy quanta. On the other hand, one also needs 
to understand the RG decoupling mentioned earlier, whereby the SU(2) theory with Np 
flavors of quark hypermultiplets flows to the SU(2) theory with Np — 1 flavors, in the limit 
that one of the quarks becomes infinitely massive. 

Both types of matching may be accomplished with the use of Weinberg's one-loop 
formula |T4] 

' l ~" ' - A(/i) . (9) 



9lM glM 

Here \i is the characteristic momentum scale of the light fields; g^ E (p) and Que (/^) are 
renormalized gauge couplings of the low- and high-energy theories, respectively; and A(/u) 
is a finite correction coming from one-loop contributions of heavy particles to the gauge 
self-energy E. Thus #he(aO is extracted from the complete set of one-loop contributions to 
E, whereas for (/ LE (/x) only light particles are permitted on the external and internal legs. 
For example, in the MS scheme, (? H e(aO and (7 LE (//) are defined in D dimensions in terms 
of the bare couplings g HEB and g hEB in the standard way: 

g HEB V 2 ~ D/2 = 9nM - b HE glM [j^i + |7b - \ lo § 47r ) , ( 10 ) 

g^ B y 2 ' ' 2 = gM - h E glM + H - I lo s 47r ) > ( n ) 



where 6 H e and 6 LE are the one-loop coefficients of the corresponding /^-functions. Note 
that when the low-energy theory is supersymmetric pure U(l) gauge theory, g^ E receives 
no perturbative corrections and is scale-independent in the limit D — > 4. 
The quantity A has the generic form 

A(m) = C 1 + C 2 log^+C 3 log^ + C4log^ , (12) 

H [i [i 

where m v , rrif and m s are the masses of the heavy vector, fermion and scalar particles 
that have been integrated out. The C, are group-theoretic constants that are tabulated 
by Hall (see Appendix 1 of ||16|| , in which an error in |TJ[] is corrected). In particular G\ 



and C2 appear only when the heavy vector particles have been integrated out; similarly, 
C3 appears when there are heavy fermions, and C4 corresponds to the heavy scalars. 

The original calculations of [|14|-|16] were performed in the dimensional regularization 
with MS scheme (DREG) . However, they can be easily translated into the supersymmetry- 



preserving dimensional reduction with MS scheme (DRED) [|T7|], or into the Pauli-Villars 
(PV) scheme. The one-loop relations between the coupling constants in different schemes 



for the SU(N C ) gauge group can be found for example in [p~8|JT 



S'pvO) #dred(» #dreg(» + 487r 2iVc ' 13 

independently of Np- It turns out that for the case of spontaneous symmetry breaking 
[SU(2), Np] — > [Z7 (1), 0] considered in detail presently, 

Cf v = Cf RED = . (14) 

In contrast, in the case of the heavy flavor decoupling [577(2), Np] — > [577(2), Np — 1], 
vector particles do not decouple and C\ = by definition. 

Now let us apply this formalism, in turn, to these two cases of interest. 

2.1. Spontaneous symmetry breaking 

We consider the Coulomb branch of iV = 2 SQCD, in which only the adjoint Higgs 
A = A a r a /2 acquires a VEV, say in the r 3 direction: (A) = ar 3 . (We have adopted 
here the VEV normalization conventions of JIJ; the translation formulae to the original 
normalization of |19[ used in our previous work are assembled in Appendix A.) The SU(2) 
component of the adjoint N = 2 supermultiplet that is parallel to the VEV then remains 
massless, while the components oc r 1 or r 2 acquire a mass M\y = Vo|a|. In addition there 
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are 2Np quark multiplets Qf and Q/, / = 1, • • • , Np, in the fundamental representation 
of the gauge group. The '1' and '2' color components of these multiplets acquire masses 
\^/2a + rrif\ and \^/2a — m,f\, respectively, as can be seen from a tree- level examination of 
the N = 2 invariant superpotential |JJ 

N F 

w = J2 >fiQt *Qf + mfQfQf . (15) 
/=i 

Here $ is the N = 1 adjoint chiral superfield whose lowest component is A; color indices 
are suppressed. 

In the PV or DRED schemes one then has 

w \ ^ M w ^ v^/, \V2a + mf\ , |\/2a — rafk ,.,„. 
A(/x) = C adj log — ^ + C f und V log J 11 + log ^) , 16 

where C adj and C fund are numerical constants. Setting fx = Myy for simplicity and extract- 
ing C fund = 1/167T 2 from Ref. fl6" 1, one finds 



1 A^F 



16tt' 

/=i 



2a 2 — m 2 



8a 2 



(17) 



so that Eq. (|9|) becomes 

1 hoop 4-JVj. , /M w \ , iV F i ^ 

In this expression <7 e g = g LE denotes the effective U(l) coupling constant, A^ F is the PV 
or equivalently the DRED dynamical scale (we drop the PV superscript henceforth) 



A^, = n b exp[-8ir/gM] = V 6 exp[-87r7^ MD (/i)] , (19) 

and we have used the fact that the (negated) coefficient of the /^-function for these models 
is b = 2N C — Np = 4 — Np. We make the following comments: 

(i) Equations (|T8"D-(|r5|) extend to Np > the case of Np = considered by Finnell 
and Pouliot, who obtained simply |L3| 



This is referred to as the absence of threshold corrections. 

(ii) As usual in a supersymmetric chiral theory l/g 2 ^ * s the imaginary part of a 
holomorphic complexified coupling constant 

87TZ 0eff 

Teff = -o- H ■ (^l) 

#eff T 

This implies that Eq. ([18[) may be analytically continued away from the imaginary axis, 
as follows: 

. N F 



r eff ^ (iM(4-N F )log(^) + ^log2 - f £ log^-m^a 2 ) . (22) 



Heavy flavor decoupling 



Next we consider the case in which a single hypermultiplet flavor becomes very heavy 
(say, m — ^ oo) and decouples from the spectrum, leaving behind N = 2 SQCD with one 
fewer flavor. Choosing (i = niN F we find that A(m7v F ) = as follows from Eq. (|T2|) with 
Ci = C*2 = 0; consequently 

9N 2 F (mN F ) = 9N 2 F -i( m N F ) (23) 



in the PV or DRED schemes. Rewriting this relation in an RG-invariant way using (|19|), 
one obtains 

m Np ■ A Np - A Np _ 1 . (24) 

The appropriate double scaling limit is therefore defined by rriN F — > oo and An f — > with 
the product on the left-hand side of (p4|) being held fixed. 

In the remainder of this section we will relate the dynamical scales A^ that appear 
in the Seiberg-Witten elliptic curves, on the one hand, to the PV or DRED scales (|19D , 
on the other hand. This is done by comparing Eqs. (|2^)-(P^) to the explicit solutions 
proposed in [Ijj. 

2.3. Relating the Seiberg-Witten scheme to the PV or DRED schemes 
The Seiberg-Witten elliptic curve for the Np = theory is |1| 

y 2 = x 2 (x-u) + ^(A S o W ) 4 x . (25) 
9 



The curve defines a dynamical scale Aq w in a particular scheme, the "Seiberg-Witten 
scheme" ; the superscript SW is introduced to distinguish it from Ao in the PV or DRED 
schemes. The correspondence between the SW and the PV schemes for Np = has 
been examined by Finnell and Pouliot |JL3|| , using Weinberg's matching formula. From the 
elliptic curve (|25|) , they calculate r e g ~ (2i/7r) log (8w/(Aq w ) 2 ) valid in the semiclassical 
regime, u ~ 2a 2 . Comparing this result with Eq. ([22]) with Np = then yields ]R 



(AD 4 = 4A*. (26) 

While this derivation was ostensibly performed at the one-loop level, it is well known that 
such relations between A's defined in different renormalization schemes are actually one- 
loop exact |[21|| . regardless of the presence of supersymmetry. (In contrast, the definition 



( |19| ) of the dynamical scale itself is one-loop exact under the RG of the Wilsonian effective 
action, but only because of supersymmetry pt)| , 



Next we consider the cases < Np < 4. Let us introduce the symmetric polynomials 
in the masses: 

N F N F N F 



The curves are then given by [plU 



Np — l 



yl p) = x 2 (x-S) + iy^(A-) 4 -^x-l(A-) 8 - 2 ^ J2 M^\x-u) N ^- s 

5=0 

(28) 

Seiberg and Witten simply equate u = u for Np = 1,2,3 but, as reviewed earlier, for 
Np = 3 this should be corrected to Eq. (|4|). It is easily checked that in the decoupling 
limit mN F — > oo, one obtains the desired result yf Np) — > 2/ ( 2 jv F -i) if an d only if one makes 
the identification 

m Np ■ (A s Z) 4 ~ Nf = • (29) 



5 A technical aside: In the absence of quark masses we have identified alternative curves for 
both Nf = 2 and Nf = 3 with the desired singularity structure. They are y 2 = x 2 (x — u) — 
^(A2 W ) 4 (x — u/9) and y 2 = x 2 (x — u) — ^-(A3 W ) 2 n 2 , respectively. Consistent with the physical 
arguments in 0, in the former case the discriminant A(u) correctly factors into a product of two 
double roots, while in the latter case A(u) is the product of a simple root with a quartic root. 
However, in each case A(n) has the wrong singularity structure once mass terms are added, which 
can be seen most easily in the case that all masses are set equal. 



10 



Notice that this is the same recursion relation as for the PV and DRED schemes, Eq. (S3). 
From the Np = "boundary condition" ( |26|) , one immediately derives the dictionary 
between schemes given in Eq. (|2|) above. 

Before passing to the case Np = 4, we note that the perturbative, one-loop structure 
of the effective U(l) coupling contained in Eq. (|22| ) has rarely appeared correctly in the 
literature. This expression has recently been confirmed by ||23|1 ; rather than invoke the 
Weinberg prescription, these authors extract the result directly from the elliptic curves, 
for arbitrary N c . 

3. Resolving the discrepancies in the Np = 4 solution 

We now turn to the interesting case Np = 4. In this model the /3-function vanishes 
and no dynamical scale is generated. Note that it is trivial to extend the RG matching 



relation (B4J) to this case. The relation (23) still holds when \i = — > oo; consequently 



m 4 exp(-8n 2 /g 2 4 ) = A 3 (30) 

using Eq. ([19|). Here, and henceforth, g^g) = #4 is the scale-independent microscopic cou- 
pling in the PV or DRED scheme!; it combines with #4 to form a single scale-independent 
PV or DRED holomorphic parameter 

8tH 64 , * 

r = — H . (31) 

Before turning to the elliptic curve, let us revisit the three discrepancies with instanton 
physics that we wish to address. The first is the discrepancy between Eqs. ([|) and (f|) 
which relate and r. (As above, we will use the superscript (0) to denote the massless 
case.) The one- loop constant cq in Eq. (ffi) was not considered in Ref. M. In light of the 



6 One may ask why, in this finite theory, it is nevertheless necessary to specify a scheme. As 
mentioned earlier, this is because the finiteness is due to an "00 minus 00" cancellation between 
divergent diagrams, which is intrinsically ill defined. For this reason we cannot simply construct 
physical quantities directly from the bare tree- level coupling gs, as would be natural to do if all 
individual graphs converged. That gB cannot be a scheme- independent physical quantity can be 
seen from Eq. (|l0]), which applies equally to the DRED and DREG schemes, and which would 
then imply for Nf = 4: <7dred(/x) = (?dreg(/u) = gs g 2 ^ ^ 2 in contradiction to Eq. (|l3|). 
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preceding discussion, we can now read off the vaiue cq = 4 log 2 from Eq. (|22|), which is 
easily extended to the case Np = 4: 

4' 4 

r eff = \?»{a) = r + -log2 - V log (l-m 2 /2a 2 ) + 0(q) . (32) 

f=i 

An independent derivation of this 1-loop relation (which tests our proposed reinterpretation 
of the massive curve) is discussed below. 

The second and third discrepancies involve properties of the massive Np = 4 curve 
(Sec. 16.3 of 0). Consider the RG decoupling property when the quark mass grows 
large. In the double scaling limit the curve indeed collapses to the Np = 3 curve (|28|); 
however, as reviewed in Appendix B, this limiting behavior requires the identification [0J 

64m 4 exp(-87r 2 /(^ w ) 2 ) = A| w = 4A 3 (33) 

where the second equality follows from Eq. (0) . This apparently contradicts Eq. ( ]3Ti| ) above; 
more precisely it means that the coupling gf w used in |jj cannot be the expected DRED 
quantity g±. And thirdly, the Np = 3 curve that one flows to in this way has u — u rather 
than the shifted definition (|j). 

As the reader can anticipate, we will posit that the parameter r that appears perva- 
sively in the Seiberg-Witten curve for Np = 4 should really be identified with the effective 
massless U(l) coupling r^, Eq. (||), rather than with the microscopic SU{2) coupling r, 
Eq. (|3lD. By definition, this reinterpretation resolves the first of the three discrepancies. 
Pleasingly it also resolves the second discrepancy: the constant factor Co = 4 log 2, when 
exponentiated, precisely compensates for the factor of 16 mismatch between Eqs. (^3|) and 
(pop. (A more stringent test of our proposal is discussed below.) Finally, the third discrep- 
ancy will be resolved by altering the relation proposed in n between u and u; however, 
there is an infinite ambiguity in this procedure that we do not know how to eliminate. 

We now review the Seiberg-Witten curve (Sec. 16 of 0), starting with the massless 

case: 

( y (o)) 2 = x 3_ i^ (Tgw)x ~2_ i^ (Tsw) ~3 = wwwwwm^ W p=x- ei (r sw )u. 

(34) 

Here (72 and (73 are rescaled Eisenstein series. The cubic roots may be defined in terms 
of ^-functions; they have the semiclassical expansion 



ei ( r ) = | + 16q 2 + 16q 4 + 0(q 6 ) , e 2 (r) = -§ - 8q - 8q 2 - 32q 3 - 8q A + 0(q 5 ) , (35) 
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with e 3 = — e\ — e 2 . As noted in |1J, the e; are not strictly speaking modular forms of 
SL(2,Z). Rather, they are weight-two modular forms of three different conjugate sub- 
groups of £X(2, Z); under the action of the full group they permute amongst themselves. 
The curve ( p^j ) is well known in the math literature JMJ . It is designed so that if the VEVs 
a and cld are extracted in the standard way as periods of the curve JIJ, 

da \[2 f dx dan V% f dx 

(36) 



du 8tt } 1x du 8tt J l2 yt°) 

then one has simply 

a = \Ju/2 , a D = r sw a . (37) 

As expected, these are the defining equations of a classical free field theory, with JF (0) (a) = 
r sw a 2 and u = 2a 2 . 

Seiberg and Witten make the two further assumptions r sw = r and u = u. The first 
of these assumptions contradicts Eq. (|J); instead, we will assume r sw = as explained 
above. As pointed out in [|J, this, in turn, invalidates the second assumption as well; 
instead, one must take 

u = u- (dr^/dr)' 1 . (38) 

This latter redefinition is specific to the massless model; it follows directly from the in- 
stanton version p5| , |T0| ,|5f of Matone's relation, which for four massless flavors reads, 



QJT(O) 

u = 27ng-^- , (39) 

combined with Eq. (|37|). 

We turn finally to the massive curve. Setting = — ej , one has |1| : 

V 2 = W 1 W 2 W 3 + e 12 e 23 e 3 i (W 1 T 1 e 23 + W 2 T 2 e 31 + W 3 T 3 e 12 ) - e\ 2 e 2 23 e 2 3l N . (40) 
Here W { = W} 0) - e 2 R, where R, N and T; are symmetric polynomials in the four masses: 
R = \H m * ' N = m 2 m 2 m 2 - ^ J] m 2 m) + ^ ^ mf , 

' i>j>k i^j i 

1 x 1 . 1 - , \ \ v^-v 

Ti = — > m 2 m 2 > mf , T 2 = I [ m, > m 2 m 2 H > mf 

1 12 24 4^ 2 11 24^ 1 3 48 ^ 1 

with T3 = — T\ — T 2 . Seiberg and Witten observe that the coefficients of this curve are 
SL(2, Z) modular invariants, provided that the aforementioned permutations amongst the 
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ej are accompanied by the same permutations acting on the Tj (this is referred to as SO (8) 
triality). They view this as strong circumstantial evidence that the dyon spectrum itself 
is SL(2, 2Hj) invariant. 

In restoring agreement with instanton calculations we need not tamper with these 
important SL(2, Z) properties. Rather, we will reinterpret the intrinsic parameters of 
the curve. In particular we will continue to take ej = ^(r^) where is given by the 
instanton series It remains only to relate the parameter u (which enters through the 
Wi) to the physical quantum modulus u = (Tr A 2 ). Dimensional analysis, 0(4) symmetry, 
smoothness in the masses, and the above-given massless limit suggest the following generic 
relation: 

c1t {0) _i 

u = u-(-p-) + R- Yl ' 9 = exp(z7rr). (42) 

(This expression is to be compared with the Seiberg-Witten proposal JIJ] 

u = u - \e x {T)R (43) 

which is already faulty in the massless limit.) The absence of odd instanton contributions 
is due to the discrete %2 symmetry discussed earlier, the mass parameter R being even 
under this symmetry. The numerical coefficients a n in (^2j) may be constrained using a 
variety of physics considerations explained in Appendix B. We find that a = —1/3, and 
«2 = 37/(3 3 2 5 ) ; the latter value disagrees again with the proposed relation (|43|), but is 
necessary to recapture the shifted definition (f|) of u in the decoupling limit. 
We conclude with the following comments: 

(i) The higher-instanton contributions a n with n = 4, 6, • • • remain completely un- 
determined, as do the constants c n with n = 4, 6, • • ■ in the relation (||) between r e ( ^ 
and r. When masses are incorporated into Eq. @, Matone's relation fl,|7|,|lOH may give 
an interesting correspondence between the two series (as it already does in the massless 
case [j5|). 

(ii) Our lack of complete knowledge of the relation between u and u does not actu- 
ally affect the low-energy effective U(l) Lagrangian. This is because a and ar> are still 
determined by Eqs. (|36| ) (with y instead of y {0) in the massive case). Both sides of these 
equalities involve the (independent) variable u rather than the unknown (dependent) vari- 
able u. The former is then eliminated in favor of a, giving a prepotential in which neither 
u nor u appears: T = T[a, u(a); {mj) = jF(a;{mi}). In other words, so long as the 
classical VEV a is considered the independent variable (and not the quantum modulus 
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u = (Tr$ 2 ) on which a depends in a presently unknown way, and vice versa), the prepo- 
tential is "known" as a function of (albeit not as a function of the microscopic r). In 
contrast, for Np < 4, both a and T are known functions of u as well. 

(iii) As a stringent test of our proposed redefinitions ei = e^r^), we have in fact 
constructed J 7 (a) as outlined in (ii), and verified that the right-hand side of Eq. (|32]) is 
indeed reproduced by the curve (paralleling Ref. [p3[ ). 

(iv) Finally we reiterate the point that the massive Np = 4 curve is still SL(2, 2Z) 
invariant, but only in terms of the rather non-intuitive quantity r^g (the effective coupling 
in the massless theory) rather than the microscopic coupling r. While the relation ([(]) 
between these parameters is currently unknown beyond the 2-instanton level, it would be 
pleasing if, in the end, the model turned out to be modular invariant in terms of r as well. 
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The work of ND was supported in part by a PPARC Advanced Research Fellowship; both 
ND and VVK were supported in part by the Nuffield Foundation; MM was supported by 
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Appendix A. Note on conventions 

Here we provide the dictionary between the normalization conventions of Ref. [ffj which 
we have adopted in this paper, and those of Ref. []TI| which we used in previous work |^,^| . 
The original unrescaled VEV definitions of []TI| are: 

&F(v) 



{A) = vr-72 , (A D ) = y D r 6 /2 , v D 



<9v 



8ivi ggg 



SQf 7T dv <9v 5 

The VEV normalizations of [PJ adopted in this paper are: 



a 



Teff 



a D 
da D 



vd 



1 dT{a) 



da 



da 

2u + ■ ■ ■ , a 



2 da 2 
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The prepotential for Np < 4 massless flavors in the two different normalizations reads: 

/ A SW \ n(4-N F ) 
71 n=2,4,6,... v " ' 

71=2,4,6,... 

Using Eq. (|2|), the relation between the instanton coefficients is determined to be: 

t(Nf) _ 9 n(6-7V F )-2 p (N F ) 
• r n — £ r n 

Appendix B. Constraints on the relation between u and u 

In this appendix we constrain some of the a priori unknown numerical constants a n 
that relate u to u as per Eq. fll!^). We use the following two considerations: 

(i) Closely following Sec. 16.3 of [|l|, we first consider the illuminating special case 
(mi, ra-2, WI3, TO4) = (0,0, m,m). Then iV = Tj = and one has simply y 2 = W1W2W3 
with Wi = x — eiU — elm 2 . As always, the critical points on the quantum moduli space are 
the values of u for which two of the three x roots coincide. In the present instance this 
means eiu + efm 2 = eju + e 2 m 2 , or equivalently since = : 

u = {eim 2 , e 2 m 2 , e 3 m 2 } ~ {|m 2 , ( — | — 8q)m 2 , ( — | + 8q)m 2 } (B.l) 

using Eq. (^5]). We have dropped terms of order q 2 as these do not survive the double 
scaling limit. Now consider the decoupling limit m — > 00. On the one hand, one expects 
a perturbative singularity at u ~ 2a 2 ~ m 2 which corresponds, physically, to two quark 
multiplets becoming massless (see Eq. (|3"2])). On the other hand, the leftover model after 
the two heavy flavors decouple is the massless Np = 2 theory, which has singularities at 
u = ±|(A| W ) 2 (see Eq. (£§)). In sum, 

u ~ {m 2 ,-i(Ar) 2 ,|(AD 2 }- (B-2) 

Equating (|B.2| ) with ( p.l| ) forces «o = —1/3; and furthermore 8q ( ^ m 2 = |(A2 W ) 2 which 
is precisely consistent with the DRED recursion relations (|29|), ( j30"D and ( |3~3"D (with g| w = 
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(ii) Returning to the case of four generic masses, we next consider the double scaling 
limit 7Ti4 — > oo, q^ — > with the product q°^ fixed at A| w /64 as per Eq. (|33|). Naively, 
we would like the curve (f|0|) to collapse to the Np = 3 curve defined by Eqs. (|28|) and (f|). 
Instead, in this limit the right-hand side of fl4C| ) diverges badly. However we can exploit 
the fact that the variable x is just a dummy of integration (see Eq. (|3~6D); an appropriate 
shift in x eliminates this divergence and guarantees a smooth RG limit. Accordingly we 



which obscures the SL(2, 2Z) properties of the curve, but preserves the Z 2 properties. 
The explicit factors in (|B.3|) have the following genesis: the constant 2/9 eliminates the 
large-mass divergence; the constant 1/3 guarantees that the cubic terms in x and u will 
have precisely the form x 2 (x — u) dictated by Eq. fl28|); and finally the factor (1/32 — a 2 /3) 
eliminates the (A^) 3 mim2m3 term that otherwise generically appears, again in order 
to harmonize with Eq. (p8|). (This term is not forbidden by any symmetry; its absence 
from the = 3 curve ( |28D is convention dependent, as is the form of the cubic term.) 
Substituting Eq. ( p.3| ) into the Np = 4 curve (|IjD and taking the double scaling limit then 
indeed reproduces the Np = 3 curve (^8|), provided that 



let 



X 



x + (l + (±-la 2 )q 2 + 0(q*))R - (\ + 0{q*))u 



(B.3) 



(l + 32a 2 )2 



-10 



= U . 



(B.4) 



From the value of uo quoted in Eq. (§j) one deduces a 2 = 37/(3 3 2 5 ) . 
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